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ABSTRACT 

We present moderate (1") and high resolution (0.2") observations of the CO 
(2-1) emission at 43 GHz, and radio continuum emission at 1.47 GHz, from 
the z = 4.7 QSO BRI 1202-0725 and the z = 4.4 QSO BRI 1335-0417 using 
the Very Large Array. The moderate resolution observations show that in both 
cases the CO emission is spatially resolved into two components separated by 
1" for 1335-0417 and 4" for 1202-0725. The high resolution observations show 
that each component has sub-structure on scales ~ 0.2" to 0.5", with intrinsic 
brightness temperatures > 20 K. The CO ladder from (2-1) up to (7-6) suggests 
a high kinetic temperature for the gas (T^in ~ 70 K), and a high column density 
(N(H 2 ) — 10^^ cm“^). In both sources the continuum-to-line ratio: — ~ 

-'"CO(l-O) 

335. All these characteristics (brightness temperature, excitation temperature, 
column density, and continuum-to-line ratio) are comparable to conditions found 
in low redshift, ultra-luminous nuclear starburst galaxies. We hnd that the CO 
emitting regions in 1202-0725 and 1335-0417 must be close to face-on in order 
to avoid having the gas mass exceed the gravitational mass, implying perhaps 
unreasonably large rotational velocities. While this problem is mitigated by 
lowering the CO luminosity-to-H 2 mass conversion factor (X), the required X 
values become comparable to, or lower than, the minimum values dictated by 
optically thin CO emission. We considered the possibility of magnihcation by 
gravitational lensing in order to reduce the molecular gas masses. 

Subject headings: radio continuum: galaxies — infrared: galaxies — galaxies: 
active, distances and redshifts, starburst, evolution, radio lines 



1. Introduction 


The recent discoveries that the large majority of spheroidal galaxies in the nearby 
universe contain massive black holes, and that the black hole mass correlates with the mass 
of the spheroid of the parent galaxy, have led to the hypothesis of co-eval formation of 
massive black holes and galaxy spheroids, perhaps occurring in merging galaxies at high 
redshift (Richstone et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 2000; Franceschini 
et al. 1999, Blain et al. 1999, Kauffmann & Haenelt 2000). We have undertaken an 
extensive study of high redshift QSOs in order to address the interesting question of co-eval 
black hole and spheroidal galaxy formation. The large samples of high ^ QSOs coming 
from wide held surveys such as the Sloan Digital Sky Survey (Fan et al. 2001), and the 
Digitized Palomar Sky Survey (DSS; Djorgovski et al. 1999) have greatly facilitated the 
study of high redshift QSOs. Our studies involve: (i) searches for thermal dust emission at 
mm-wavelengths from a large sample oi z > 3.7 QSOs (Carilli et al. 2001a; Omont et al. 
2001), (ii) high resolution imaging at cm-wavelengths of the non-thermal radio continuum 
emission from these sources (Carilli et al. 2001a, Carilli et al. 2001b), and (hi) observations 
of the CO line emission from selected sources with large infrared luminosities, as inferred 
from the mm-continuum observations (Carilli, Menten, & Yun 1999). 

The interesting result from the mm-continuum surveys of high redshift QSOs is that 
30% of the sources are detected in surveys with flux density limits of 1 to 2 mJy at 250 
GHz (Carilli et al. 2001a; Omont et al. 2001), corresponding to infrared luminosities 
> 10^^ Lq and dust masses > 10® Mq. Radio continuum studies of these sources show that 
the mm-to-cm spectral energy distributions (SEDs) of most of the sources are consistent 
with the radio-to-far IR correlation found for nearby star forming galaxies (Condon 1992; 
Carilli 2001a, b; Yun et al. 2000). If the dust is heated by star formation, the implied star 
formation rates are of order 10® Mq year“^. On the other hand, the global SEDs for these 
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sources from cm-to-optical wavelengths are not out of the range dehned by lower redshift, 
lower luminosity QSOs (Carilli et ah 2001a,b; Sanders et ah 1989), and hence the case for 
dust-heating by star formation in these sources is by no means secure. 

Searches for CO line emission from dust emitting QSOs at 2 : > 4 have resulted in 
the detection of molecular line emission from four sources to date: 1202-0725, 1335-0415, 
0952-0115, and 2322-1-1944 (Ohta et ah 1996; Omont et ah 1996b; Guilloteau et al. 1997, 
1999; Cox et al. 2001). The implied gas masses are large (> 10^*^ Mq), leading some to 
speculate that star formation may be inevitable (Omont et al. 2001). Moreover, in at 
least one case (1202-0725) the molecular line and mm continnnm emission is known to be 
spatially extended, with a bright emission region well separated from the optical QSO. 

For sources at 2 ; > 3.9 the CO(2-l) transition is redshifted into the 43 GHz band of 
the Very Large Array (VLA), thereby allowing for sensitive observations to be made of the 
low order transitions from these sources. The VLA also allows for sub-arcsecond angular 
resolution observations to be made, which are critical to the physical interpretation of the 
systems. 

In this paper we present high resolution observations of the CO(2-l) emission from 
1202-0725 at z = 4.7 and from 1335-0415 at z = 4.4, along with radio continuum 
observations at 1.4 GHz of 1202-0725. These data are combined with observations of higher 
order CO transitions in an attempt to understand the physical conditions in the molecnlar 
line emitting regions of high 2 ; QSOs. We use Ho = 65 km s“^ Mpc“^, Qm = 0.3 and 
Oa = 0.7. 
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2. The sources 

By selecting for very red, point-like optical sources, McMahon (1991) has identihed a 
large sample of z > 4 QSOs from the Automatic Plate Measuring survey (Irwin, McMahon, 
& Hazard 1991), including the two sources BRI 1202-0725 at z = 4.7 and BRI 1335-0415 
at z = 4.4. Optical spectra of both of these sources show strong associated Lya absorption 
(Storrie-Lombardi et al. 1996). 


2.1. 1202-0725 

Omont et al. (1996) and Ohta et al. (1996) detect high order CO line emission 
from 1202-0725, as well as thermal continuum emission from warm dust at 1.35mm. The 
continuum and line emission consist of a double source with an angular separation of 4". 
The total flux density of the source at 1.2mm is 12 ± 3 mJy, implying a far IR luminosity 
of Lfir = 4.2 X 10^^ Lq, assuming a spectral energy distribution (SED) typical for an 
ultra-luminous infrared galaxy (ULIRG; Tfir ~ 10^^), where Tfir is dehned as in Helou et 
al. (1988) and Condon (1992), ie. the integrated luminosity between rest frame wavelengths 
of 42 to 122 /rm. The southern source in 1202-0725 comprises about 65% of the total. 
Observations of this source with ISO at mid- to far-IR wavelengths constrain the dust 
temperature to be ~ 70 K for a dust emissivity index /? = 1.5 (Leech, Metcalfe, & Altieri 
2001 ). 

The CO observations of Omont et al. (1996a) and Guelin et al. (2001) indicate a 
difference in the CO line widths for the northern and southern components. Gaussian 
fitting to the CO(5-4) line emission prohles by Omont et al. (1996) results in z = 4.6947, 
line Full Width at Half Maximum (FWHM) = 190 km s“^, and an integrated line flux 
density = 1.1±0.2 Jy km s“^ for the southern source, and z = 4.6916, FWHM = 350 km 
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s“^, and an integrated line flux density = 1.3 ± 0.3 Jy km for the northern source. 
This difference is important, since it argues against gravitational lensing as the origin for 
the observed double source and large apparent luminosity. However, Guilloteau (2001) has 
recently called the difference in CO profiles into question due to the possible contamination 
of the line emission from the northern source by continuum emission at 3mm. The data 
presented herein support the original conclusions of Omont et al. (1996a) and Guelin et al. 
(2001) for different line profiles for the northern and southern sources in 1202-0725. 

BRI1202-0725 is a radio continuum source with a total flux density of 315 ± 80 ^uJy at 
1.4 GHz (see section 4.1), and 141T15 ^uJy at 4.9 GHz (Yun et al. 2000). Extrapolating 
the centimeter spectrum to millimeter wavelengths assuming a continuous powerlaw 
spectrum implies a non-thermal contribution to the integrated emission at 40 GHz of 
80/iJy. Conversely, extrapolating the (sub)mm dust emission spectrum downward to 40 
GHz implies a thermal dust contribution of 25^Jy, assuming j3 = l.b. 

Optical imaging of 1202-0725 shows a point-like QSO with = -28.5, with faint 
near-IR and Lya emission extending about 2.4" north of the QSO (Hu, McMahon, & Egami 
1996). The optical position of the QSO given by Hu et al. is 0.5" north of the position of 
the southern mm/cm/CO component. In the following we assume that this offset indicates 
the accuracy of the relative astrometry of the radio and optical images. 

The (3cr) limit to the Ly a flux from the northern mm/cm/CO component in 1202-0725 
is ~ 1 X 10“^^ erg cm“^ s“^ (Hu et al. 196), corresponding to a limit to the Ly a luminosity 
of 2 X 10^^ erg s For comparison, the Lya luminosities of UV dropout galaxies at z ^ 3.09 
are typically 2 to 4 x 10^^ erg s“^ (Steidel et al. 2000). The line emitting object detected 
2.4" northwest of the QSO has a Lya luminosity of 4 x 10^^ erg s“^ (Fontana et al. 1998; 
Petitjean et al. 1996). If this line emission is powered by star formation, then the star 
formation rate is > 20 Mq year“^, depending on the amount of dust extinction (Fontana et 
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al. 2000). It is also possible that this line emission is powered by UV radiation from the 
QSO itself (Petitjean et ah 1996). The lack of CO and thermal dust emission from this 
position argues for the latter, while the lack of N V and C IV emission lines argues for the 
former (Fontana et al. 1998; Ohta et al. 2000). 


2.2. 1335-0417 

Guilloteau et al. (1997) have detected CO (5-4) line emission from BRI 1335-0417 at 
2 ; = 4.4074, with a line FWHM = 420 km s“^ and an integrated line intensity of 2.8±0.3 
Jy km s“^. The flux density of the dust continuum emission measured at 1.35mm with the 
Plateau de Bure Interferometer (PdBI) at 2" resolution is 5.6±1.1 mJy, while that measured 
at 1.25mm with the IRAM 30m telescope (10.6" resolution) is 10.3 ± 1.4 mJy. The implied 
far IR luminosity based on the 1.25 mm measurement is 3.1 x 1013 Lq. The 1.35mm PdBI 
observations shows marginal evidence for an extended source, with a formal Gaussian size 
of 1.0 ± 0.4", with major axis oriented roughly north-south. 

No high resolution optical images have been published of 1335-0417. The DSS shows 
an unresolved source with = -27.3, although with a pixel scale of just 1" the details 
of the source structure remain unknown. BRI 1335—0417 has been detected in the radio 
continuum, with an integrated flux density of 220 ± 43 /iJy at 1.4 GHz and 76 ± 11 /xJy at 
4.9 GHz (Carilli et al. 1999). 


3. Observations 

Table 1 lists the observing parameters for BRI 1202-0725 and BRI 1335—0417. We 


also present a re-analysis of the D conhguration observations from Carilli et al. (1999). 
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Standard amplitude and phase calibration were applied, correcting for atmospheric 
opacity at high frequency, and the absolute flux density scale was set by observing 3C 286. 
Fast switching phase calibration was employed for the high frequency observations (Carilli 
& Holdaway 1999). The 43 GHz observations were dynamically scheduled, and took place 
at night under excellent weather conditions. RMS phase variations after calibration were < 
20°. The phase coherence was checked by imaging a calibrator with the same calibration 
cycle as that used for the target sources. At all times the coherence was found to be higher 
than 85%. 

A severe limitation at the VLA for observing broad lines at high frequencies is the 
maximum correlator bandwidth of 50 MHz, and the limited number of spectral channels (7) 
when using this bandwidth and dual polarization. The bandwidth of 50 MHz corresponds 
to a velocity coverage of only 350 km s“^ at 43 GHz. These correlator limitations preclude a 
meaningful determination of the line prohle, so we chose to observe in continuum mode with 
two IFs (Intermediate Frequency) with two polarizations for each IF in order to maximize 
sensitivity to the integrated line emission (Garilli et al. 1999). Our analysis will necessarily 
assume the line widths as given by the higher order transitions. 

For the GO line observations of 1202-0725 IFl was centered on the emission line from 
the southern source, while IF2 was off-set from this by 50 MHz. Based on the (admittedly 
noisy) spectra in Omont et al. (1996a), IFl should contain most of the emission from 
the southern component, while the sum of the IFs cover the emission from the northern 
component. The continuum for 1202-0725 was investigated with an observation at 43 GHz. 
For 1335-0417 IFl was centered on the emission line, while IF2 was centered 1400 km 
s“^ off the line. Taking into account the limited width of the continuum band employed, 
we estimate that we are missing about 30% of the velocity-integrated line emission in 
1335-0417. 
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4. Results 

4.1. BRI 1202-0725 

The VLA D array observations of the CO(2-l) emission from 1202-0725 are shown 
in Figure 1, with a spatial resolution of about 2". Figure la shows the summed emission 
in IF’s 1 and 2 for the 40 GHz observations. Shown in greyscale on Figure la is the Ly 
a image of 1202-0725 from Hu et al. (1996). Again, we have aligned the position of the 
optical QSO with the peak in the radio continuum, CO, and thermal dust emission for 
the southern source. Figures lb and Ic show IFs 1 and 2, respectively, for the 40 GHz 
observations. Figure Id shows the continuum image at 43 GHz. 

Table 2 lists the results for the two components in 1202-0725. Column 3 lists the 
position of the measured CO components. Column 4 lists the 7mm continuum flux densities 
at these positions, column 5 lists the observed CO flux densities, and column 6 lists the 
velocity integrated CO emission. Column 7 lists the 1.4 GHz continuum flux densities, 
while column 8 lists the CO luminosities (in K km s“^ pc^). 

The southern CO component is clearly detected in IFl of the 40 GHz observations 
(Fig. lb), but is not seen in IF2, nor in the 43 GHz continuum image (Figs. Ic and Id). 
The northern component shows up in both IF channels for the 40 GHz observations (Figs 
la,b,c), but is not seen in the 43 GHz continuum image (Fig Id). From the 43 GHz image 
we set a 2a upper limit of 0.16 mJy for the 43 GHz continuum emission from both the 
northern and southern components in 1202-0725. 

The lack of continuum emission at 43 GHz from the northern component in 1202-0725 
is signihcant in regard to the recent possible detection of 0.8 ± 0.15 mJy of continuum 
emission at 3mm by Guilloteau (2001). The implied spectrum from 3mm to 7mm must be 
rising with a powerlaw index > 1.9. It also implies that the observed emission in IFs 1 and 
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2 at 40 GHz for the northern source is CO(2-l) line emission (Figs. la,b,c), and hence that 
the velocity prohles for the CO lines from the northern and southern sources are different. 

Figure 2 shows the results from the B array observations with a resolution of about 
0.25". The southern source in IFl (Fig 2b) appears as two unresolved, roughly equal 
components (flux densities ~ 0.41 ± 0.12 mJy), separated by 0.3". The implied lower limit 
to the (redshift corrected) brightness temperature for these components is about 25 K. 
The northern source in IFs 1 and 2 (Fig 2a) is marginally detected with a peak surface 
brightness of 0.25 ± 0.08 mJy beam“^. This implies two, or more, compact components 
with flux densities below our detection threshold, or diffuse emission on a scale > 0.5". 

Figure 3 shows the radio continuum image of 1202-0725 at 1.4 GHz, with a resolution 
of about 2". Both components are detected, with flux densities as listed in Table 2. The 
total 1.4 GHz flux density of 1202-0725 has been measured three times over the course of 
three years with values of 240 ± 40 yuJy (Yun et ah 2000), 305 ± 60 /xJy, and 390 ± 40 
/iJy. In our analysis we adopt the mean value of 315 /iJy. It is possible that the source is 
variable, although making accurate flux density measurements at this level is difficult due 
to confusion problems arising in wide-held imaging at 1.4 GHz. 


4.2. BRI 1335-0417 

The low resolution image of the CO(2-l) emission from 1335-0417, as reproduced from 
Carilli et ah (1999), is shown in Figure 4a, along with the off-line channel in Figure 4b. For 
this image the UV-data were naturally weighted in order to maximize sensitivity, at the 
expense of resolution. The rms noise on the image is 0.11 mJy beam“^ and the resolution 
is about 1.6". The CO emission appears extended north-south in this image by about 1". 


In order to investigate this extension in more detail, we re-imaged the data using 
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uniform weighting of the visibilities, which optimizes resolution at the expense of sensitivity. 
The result is shown in Figure 4c, with a rms noise of 0.14 mJy and resolution of about 1.3". 
The source is comprised of two components separated by 1.3". Gaussian htting to each 
component shows that they are unresolved, with upper limits to their sizes of about 1.1". 
Table 2 lists the component positions and flux densities. 

The position of the 1.35mm continuum peak is located within 0.2" of the southern CO 
component, and the mm continuum source also shows marginal evidence for a north-south 
extension on the scale of 1" (Guilloteau et al. 1997). The optical QSO position is within 
0.1" of the southern source position. 

Figure 5 shows the high resolution (B array) image of the CO(2-l) emission from 
1335-0417. The naturally weighted beam is roughly circular with FWHM = 0.17". 
Nothing is detected in this high resolution image to a 2a surface brightness limit of 0.18 
mJy beam“^. Non-detection at high resolution could mean that the emission is diffuse 
on scales larger than the resolution with a (redshift corrected) brightness temperature 
< 22 K. On the other hand, a number of compact components distributed over 1" with 
higher brightness temperature is certainly not precluded, eg. the data allow for four small 
(< 0.17") components each with brightness temperatures > 22 K. 


5. Analysis 
5.1. Masses 

The H 2 gas masses can be calculated from the values of L' in Table 1 assuming a value 
of X = the H 2 mass-to-GO(l-O) luminosity conversion factor in Mq (K km s“^ pc^)“^. 

A value of X ~ 4.6 is applicable to Galactic Giant Molecular Glouds (Dame et al. 1987; 
Strong et al. 1988; Bronfman et al. 1988), while for ULIRGs Downes and Solomon (1998) 
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find a value of X ~ 0.8. This calculation also requires an extrapolation from the CO(2-l) 
measurements to the CO (1-0) luminosity. Assuming constant brightness temperature and 
using a value of X appropriate to ULIRGs leads to molecular gas masses of about 5 x 10^° 
Mq for the CO emitting components in 1202-0725 and 1335-0417. 

We can also calculate the gravitational masses using the observed sizes and line widths 
for the sources, and assuming a disk of gas in Keplerian rotation. For 1202-0725 south we 
use a line width of 190 km s“^ and a radius corresponding to half the separation of the two 
components = 0.15" = 1.1 kpc, leading to an enclosed mass of 0.23 x 10^'^ sin“^(f) Mq. 
For 1335-0417 we do not have spatially resolved spectroscopy, so we adopt a line width of 
420 km s“^ and a radius of half the separation of the two components = 0.65" = 4.7 kpc, 
leading to an enclosed mass of 4.8 x lO^'^ sin“^(i) M©. 

An upper limit to i can be derived by assuming that the molecular gas mass dominates 
the total enclosed mass. Using X = 0.8 for 1202-0725 south we hnd i < 12°, while for 
1335-0417 we hnd i < 43°. These small angles lead to a signihcant problem in that the 
implied rotational velocities are then extremely large, > 400 km s“^. To obtain rotational 
velocities of ~ 250 km s“^, more typical of large spiral galaxies, would require a further 
reduction in the conversion factor to X < 0.2. However, such a low X value would violate 
the minimum mass conditions dictated by optically thin CO emission, for which X ~ 0.35 
for an excitation temperature of 50 K and assuming a Galactic CO abundance (Solomon 
et al. 1997). It may be that these systems are extremely massive, or perhaps that the 
apparent CO luminosities are magnihed by gravitational lensing (see section 6). 
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5.2. Continuum-to-Line Ratios 

We now consider the continuum-to-line ratio: —. A non-linear relationship 

''^CO(l-O) 

between this ratio and Lfir has been found by Solomon et al. (1997), in which values of 5 
to 50 have been found for Galactic Giant Molecular Glouds and for nearby galaxies with 
-hpiR < 10^°, while values of 80 to 250 have been found for ULIRGs. 


For the QSOs discussed herein we derive the L'(GO(1-0)) from L'(GO(2-l)) assuming 
the same brightness temperature for the two transitions. Gonsidering the integrated 
properties of 1202-0725 leads to: = i 2 xio^iK^km^s-i pG ~ 1335-0417 we 


hnd: 


^FIR _ 

r / 

^CO(l-O) 


9.6xlO^^K km pc^ 


^CO(i-O) 

= 323. 


5.3. CO excitation conditions 

Figure 6 shows the GO ladder for the various transitions observed in 1202-0725 and 
1335-0417. The data have all been normalized to the velocity integrated line flux for 
GO (5-4). For comparison, we have included the GO ladder observed for the best studied 
nuclear starburst M82 (Glisten et al. 1993; Mao et al. 2000), and the GO ladder for the 
integrated emission from the the Milky Way disk inside the solar radius (excluding the 
Galactic center) as seen by GOBE (Fixsen, Bennett, & Mather 1999). The excitation 
conditions for the two QSOs follow those seen in M82, with roughly constant brightness 
temperature (intensity oc frequency^) up to GO(5-4), and then roll-off to higher order 
transitions. This behavior is very different than that seen for the disk of the Milky Way, for 
which the integrated line flux peaks at the GO(3-2) line. 

We have used a standard one component LVG model simulating a spherical cloud to 
interpret the observed line ratios (see also Ohta et al. 1998). Input parameters are the 
kinetic temperature, the H 2 density, and the GO column density, the latter in the form 
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N(CO) = 3.08 X 10^® X iico x [ grtdv ] with grad V denoting the velocity gradient in 

km pc“^, n(CO) being the nominal CO number density in cm“^, and AV describing the 
FWHM in km s“^. Given the inhomogeneity of the interstellar medium in the Milky Way 
and nearby galaxies, our assumption of uniform physical conditions is crude but appropriate 
for a source with four detected lines from the main CO species and no information on rare 
isotopomers. The cosmic background was assumed to be at a temperature of 15 K. 

Due to the lack of constraints, there is a clear degeneracy in LVG modeling between 
various parameters such as the temperature and density (Giisten et ah 1993). Our 
modeling is not meant as an exhaustive analysis, but is merely representative of the types 
of conditions that can give rise to the observed line ratios. Adopting a kinetic temperature 
of 70 K, corresponding to the temperature derived from the spectrum of thermal dust 
emission (Leech et al. 2001), we hnd that the data can be reasonably £t with a model in 
which: N(CO) = 5 x 10^® cm“^, n(H 2 ) = 2 x 10^ cm“^, and grad V = 1 km s“^ pc“h 
Raising the temperature to 100 K decreases the required density by a factor two. A 
more thorough analysis of the excitation conditions in 1202-0725, including new, spatially 
resolving observations of the high order transitions will be given elsewhere (Kohno et al. 
2001, in preparation). 

In the likely case that photon dominated regions (PDR) provide a more realistic 
approach, the density must be larger than the ‘critical’ density to efficiently excite the mid- 
and high-J CO lines and to avoid subthermal excitation (e.g. Koester et al. 1994). With 
line temperatures being highest for the CO J = 2-1 or 3-2 transitions, critical densities, at 
which collisional deexcitation rates exceed radiative decay rates, are 10^ to 10^ cm“^. 
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6. Discussion 

We briefly compare the properties of 1202-0725 and 1335-0417 to those found in 
nuclear starburst galaxies seen at lower redshift. Yun et al. (2000) found that the ratio of 
radio-to-far IR luminosity for both these sources is within the range dehned for active star 
forming galaxies based on the tight radio-to-far IR correlation (Condon 1992), although 
in both cases the ratio falls at the high end of the normal range, suggesting a possible 
contribution to the radio emission from the AGN. For 1202-0725 there is marginal evidence 
that the radio emission is time variable. Variability would rule-out a starburst origin for the 
radio emission. Further radio continuum monitoring is in progress to test this interesting 
possibility. 

The measured (redshift corrected) brightness temperature for the CO(2-l) emission 
from the compact components in the southern source in 1202-0725 is > 25 K. This limit is 
consistent with values seen for the CO(2-l) emission from starburst nuclei of ULIRGs, for 
which brightness temperatures of 30 K to 60 K have been measured (Downes and Solomon 
1998). 

The excitation conditions for the CO for these two sources follow roughly those seen 
for the nuclear starburst galaxy M82, but are very different than those expected for the disk 
of a normal spiral galaxy. Assuming a Galactic abundance for CO of = 5 x 10“^ (eg. 
Wilson et al. 1986) implies an H 2 column density of order 10^^ cm“^. This value is similar 
to the molecular gas column density seen toward ULIRGs (Downes & Solomon 1998), but 
is an order of magnitude larger than that observed in the lower luminosity nuclear starburst 
galaxies M82 and NGC 253 (Harrison, Henkel, & Russell 1999, Mao et al. 2000). A low 
metallicity would suppress CO emission and would yield even higher H 2 column densities. 
It would be very interesting to measure the CO(l-O) line of the source. A direct comparison 
of CO(l-O) and 2-1 line intensities would provide signihcant information on optical depths. 
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since the 1-0 line may be optically thin, while the 2-1 line is likely optically thick. 

We have found that the continuum-to-line ratio, , is about 335 for both sources. 

-^CO(l-O) 

This ratio is at the high end of those seen for ULIRGs, and suggests a continuation of the 
non-linear trend for increasing continuum-to-line ratios with increasing far IR luminosity. 
Considering Lpm fo be a measure of star formation rate, and L'(CO(1-0)) to be a measure 
of molecular gas mass, it has been suggested that this non-linear relation might imply 
a higher star formation efficiency (= star Rate ^ higher luminosity galaxies, in 

particular for galaxies with Lfir > 10^^ Lq (Solomon et al. 1997). For dense nuclear 
starbursts a number of groups (Solomon et al. 1997; Mao et al. 2001; Weiss et al. 2001) 
have shown that the densities are such that the entire interstellar medium in the starburst 
regions may be molecular, and that the CO (1-0) emission may be dominated by this 
molecular inter-cloud medium, as opposed to being from the denser star forming clouds 
themselves. This phenomenon would contribute to the non-linear relationship between Lpm 
and L'(CO(1-0)). For 1202-0725 and 1335-0417 there is also the obvious possibility of dust 
heating by the AGN, in which case the continuum-to-line ratio cannot be interpreted in the 
context of star formation. 

Downes and Solomon (1998) show that for the nuclear starburst regions in ULIRGs 
the H 2 mass-to-CO(l-O) luminosity conversion factor, X, is a factor four or so below the 
Galactic disk value. Even for this low value of X, we hnd that the CO emitting regions in 
1202-0725 and 1335-0417 must be close to face-on in order to avoid having the gas mass 
exceed the gravitational mass, implying perhaps unreasonably large rotational velocities. 
While this problem is mitigated somewhat by lowering X even further, the required X 
values become comparable to, or lower than, the minimum values dictated by optically thin 
CO emission (Solomon et al. 1997). 


One way of circumventing this mass problem would be to assume that the source is 
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magnified by strong gravitational lensing. Magnification by a factor three or so would avoid 
unphysical values of X. The very large IR and CO luminosities, and the double structure of 
the sources in the CO line emission, and in the non-thermal radio continuum and thermal 
mm continuum for 1202-0725, could possibly indicate gravitational lensing. On the other 
hand, neither source is double at optical wavelengths, arguing against lensing, although the 
possibility of differential obscuration along the two lines-of-sight complicates this conclusion 
(Hu et ah 1996). More telling is the difference in CO(2-l) line prohles, which is difficult, 
although perhaps not impossible, to explain in the context of gravitational lensing. Also, 
double sources arise from strong gravitational lensing of very compact emitting regions 
(sizes < few pc). For extended emitting regions (> 100 pc), such as must be the case 
for the thermal CO and mm continuum emission, strong lensing will only occur if the 
extended emission regions cross a caustic in the source plane (Blandford & Narayan 1992). 
Such a phenomenon usually leads to more complex geometries, like arcs or rings, as in 
APM 0827+525 and IRAS 10214 (Lewis et al. 2001; Scoville et al. 1995). More sensitive, 
high resolution imaging at cm and mm wavelengths is required to address this interesting 
question. 

Overall, the physical conditions in the molecular gas and dust in these systems are 
similar to those observed in nuclear starbursts at low redshift, including: (i) the radio-to-far 
IR luminosity ratio, (ii) the CO brightness temperature, (iii) the CO excitation conditions, 
(iv) the CO column densities, and (v) the CO line-to-dust continuum ratio. 

The National Radio Astronomy Observatory (NRAO) is operated by Associated 
Universities, Inc. under a cooperative agreement with the National Science Foundation. 
We thank E. Hu for allowing us to reproduce the Lya image, S. Myers for discussions 
concerning gravitational lensing, and the referee for many important comments. 
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Figure Captions 

Fig. la.— The contours show the VLA image of C0(2-l) emission from 1202-0725 at 
= 4.691. The spatial resolution (FWHM) is 2.6" x 1.8" with the major axis position angle 
= 0°, and the rms noise on the image is a = 0.07 mJy beam“^. This frame shows the sum 
of the two IFs with a mean frequency of 40.510 GHz and a total bandwidth of 100 MHz. 
The contour levels are: -0.13, 0.13, 0.26, 0.39 mJy beam“^. The greyscale shows the narrow 
band image of the Lya emission (Hu et ah 1996). The crosses in this and subsequent 
images of 1202-0725, mark the two peak positions of the CO emission in this image. 

lb. — The image of IFl (40.485 GHz) and a = 0.1 mJy beam“^. The contour levels are: 
-0.19, 0.19, 0.38, 0.57, 0.76 mJy beam“h 

l c. — The image of IF2 (40.535 GHz) with the same contours as lb and a = 0.1 mJy 
beam“^. 

l d. — The 43 GHz continuum image of 1202-0725 with a resolution of 3" x 1.6", with the 
major axis position angle = —20", and a = 0.08 mJy. The contour levels are: -0.2, -0.1, 0.1, 
0.2, 0.3 mJy beam“^. 

Fig. 2a.— The VLA image of the CO(2-l) emission from the northern component in 
1202-0725 using the sum of IFs 1 and 2 at 40 GHz at a resolution of 0.30" x 0.21" with the 
major axis position angle = 11", and a = 0.08 mJy beam“^. The contour levels are: -0.2, 
-0.1, 0.1, 0.2, 0.3 mJy beam“L 

2b.— The VLA image of the CO(2-l) emission from the southern component in 1202-0725 
for IFl at 40.485 GHz, and a = 0.12 mJy beam“^. The contour levels and resolution are 
the same as Figure 2a. 


Fig. 3.— An image of the radio continuum emission from 1202-0725 at 1.4 GHz made 
with the VLA at a resolution of 2.0" x 1.5", major axis position angle = -30°, and cr = 16 
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yuJy beam The contour levels are a geometric progression in the square root of two with 
the hrst level being 35/iJy beam“^. 

Fig. 4a.— The contours show the VLA image of CO(2-l) emission from 1335-0417 at 
2 : = 4.4074. The spatial resolution is 1.4" x 1.3" with major axis PA = 39°, and a = 0.11 
mJy beam“^. This frame shows the first IF at 42.635 GHz and a bandwidth of 50 MHz. 
The contour levels are: -0.24, -0.12, 0.12, 0.24, 0.36, 0.48, 0.50 mJy beam“^. The crosses in 
this and subsequent images, mark the positions of the two peaks as seen in Figure 4c. 

4b. — The same as 4a but for the second IF at 42.835 GHz. 

4c. — The same as 4a but using uniform weighting of the visibilities, leading to a spatial 
resolution of 1.4" x 1.3" with major axis PA = 39°, and a = 0.14 mJy beam“^. The contour 
levels are: -0.26, -0.13, 0.13, 0.26, 0.39, 0.52 mJy beam“^. 

Fig. 5. — An image of the GO(2-l) emission from 1335-0417 at a spatial resolution 
of 0.17", and a = 0.09 mJy beam“^. This frame shows the hrst IF at 42.635 GHz and a 
bandwidth of 50 MHz. The contour levels are: -0.28, -0.14, 0.14, 0.28, 0.42 mJy beam“F 

Fig. 6. — The GO ladder for the two high redshift QSOs in this paper and using data 
from Omont et al. (1996), Ohta et ah (1996), and Guilloteau et ah (1997) for the other 
transitions. The open squares are the results for BRI 1202-0725. The open triangles are 
those for 1335-0417. The solid triangles are the data for the GO ladder for the integrated 
emission from the the Milky Way disk inside the solar radius (excluding the Galactic center) 
as seen by GOBE (Fixsen, Bennett, & Mather 1999). The open circles are the results for 
the starburst nucleus of M82 (Glisten et al. 1993; Mao et al. 2000). The ordinate is the 
velocity integrated line flux density, normalized to the GO(5-4) line, except for the Milky 
Way, for which the values are normalized at GO(3-2). The long dash line shows the line 
strengths increasing as frequency squared. The short dash line shows an LVG model with 
Tkin = 100 K and n(H 2 ) = 1 x 10^ cm“^, and the dotted line shows the LVG model with 
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Tkin = 70 K and n(H 2 ) = 2 x 10^ cm 



Table 1. VLA Observations 


Source 

Date 

Configuration 

Frequency 

GHz 

1202-0725 

March 1999 

D (1km) 

40.485, 

40.535 

1202-0725 

May 2001 

B (10km) 

40.485, 

40.535 

1202-0725 

November 2001 

D 

43.315, 

43.365 

1202-0725 

August 1999 

A (30km) 

1.365, 

1.415 

1335-0417 

December 1999 

B 

42.635, 

42.835 

1335-0417 

March 1999 

D 

42.635, 

42.835 
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Table 2. Observed Parameters 


Source 

z 

Position 

J2000 

S(7mm) 

mjy 

S(CO(2-l)) 

mJy 

S AV 

Jy km s~^ 

S(1.4 GHz) 

mJy 

L’(CO(2-l)) 

Xl0^° K km pc^ 

1202-0725 South 

4.695 

12 05 23.12 -07 42 32.9±0.3" 

0.16 ± 0.08 

0.77 ± 0.10 

0.23 ± 0.04 

0.16 ± 0.016 

5.5 ± 1.0 

1202-0725 North 

4.692 

12 05 22.98 -07 42 29.9±0.3" 

-0.13 ± 0.08 

0.44 ± 0.07 

0.26 ± 0.05 

0.23 ± 0.016 

6.2 ± 1.3 

1335-0417 South 

4.407 

13 38 03.40 -04 32 35.4±0.3^^ 

-0.03 ± 0.12 

0.67 ± 0.14 

0.26 ± 0.06 

0.22 ± 0.043 

5.7 ± 1.1 

1335-0417 North 

4.407 

13 38 03.42 -04 32 34.1±0.3" 

0.07 ± 0.12 

0.45 ± 0.14 

0.18 ± 0.06 

- 

3.9 ± 1.1 
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